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We studied the phase progression properties of normal modes in a 
conical waveguide in order to develop techniques for analysis of nndtimode 
microwave antennas. We found that the large-order asymptotic expansions 
of Bessel functions developed by F. W. J. Olver are most appropriate 
for such calcidations by virtue of their simplicity and uniformity with 
respect to argument. These expansions are applied to a?ialysis of the conical 
TE n and conical TM U modes and, in addition, to an examination of 
the " quasi-cylindrical approximation" in which the conical waveguide 
is regarded as a cylindrical waveguide with gradually changing cross section. 

I. INTRODUCTION 

For most applications to microwave communication systems, wave- 
guides are designed in such a way that only the dominant mode can 
propagate. This has been the case principally for practical reasons, as 
evidenced for example by problems encountered in the development 
of millimeter-wave systems using the higher-order TE^ mode. 1 Since 
the waveguide in this case must be oversized, small geometrical asym- 
metries due to errors in fabrication, bends, and other structural perturba- 
tions cause coupling to unwanted modes, which can result in a sig- 
nificant degradation in performance. 

On the other hand, there has been considerable interest during the 
past several years in techniques winch require the controlled excitation 
of higher-order modes combined with the dominant mode in, for ex- 
ample, a conical waveguide. Two such applications have been the 
TE< - TMox* precision autotrack system for the Telstar® satellite, 2 
and the TE^ — TM^ dual mode conical horn 3 which has been suggested 
as a primary feed configuration for low-noise satellite communication 
antennas. 



* The notation TE< n or TM£ n will be used to designate conical waveguide 
modes. 
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A common feature of these techniques is the necessity for maintain- 
ing a high degree of phase coherence among the various modes of 
propagation. It is therefore required, in order to design such systems 
and predict the effects of frequency, temperature and structural varia- 
tions, to accurately determine the phase progression properties of the 
guided wave fields. 

The paper's contents may be summarized as follows: We first de- 
scribe the conical waveguide modes, which are vector eigenfunctions 
of Maxwell's equations in what is essentially spherical geometry. We 
note that although these solutions are well known in principle, the 
actual computation of their phase progression properties is not 
straightforward. It is, therefore, necessary to consider the problem of 
numerically evaluating both the eigenvalues and vector eigenfunc- 
tions so that we can apply these results to actual antenna problems. 
In order to do this we utilize certain uniform asymptotic expansions 
due to F. W. J. Olver 4 which are found to be well suited to such cal- 
culations. We thereby observe that a very common method of deter- 
mining phase progression, which might be termed the quasi-cylindrical 
approximation, is not particularly accurate, and the errors associated 
with this method are evaluated. 

Rationalized MKS units and the (suppressed) harmonic time de- 
pendence e~ /u)f will be used throughout. 

II. MODES IN A CONICAL WAVEGUIDE 

The normal modes characteristic of a conical waveguide are derived 
in the usual manner by finding separable solutions to Maxwell's equa- 
tions in spherical coordinates, subject to the boundary condition that 
the components of electric field tangent to the lateral surface must 
vanish. The solutions thus derived may be partitioned into two types, 
TE* modes for which the electric field is transverse to the direction 
of propagation (the r-direction), and TM < modes for which the mag- 
netic field is transverse to the direction of propagation. In terms of 
the coordinate system shown in Fig. 1, the components of electric field, 
for example, are given by, 5 

TM< 

\jcr) 

£ »™ - idk t ^ H ^ kr)] To p "- (cos e)e '" (1) 
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E™ = 



Fig. 1 — Conical horn geometry. 



TE* 



E™ = 

e: e = 



sin By/kr 

iB 



Hl\\{kr)P m Xcos 6)e" np 



(2) 



-.HlV i (kr)^[P m Xco B B)]e in r 



In (1) and (2), P m „ (cos 6) denotes the associated Legendre function 
(w = integer) and H[ l) (x) represents the Hankel function of the first 
kind, corresponding to outgoing waves under the assumed time de- 
pendence e~' wt . The constant k is the free-space wave number. 

The eigenvalues v and /x are found as solutions of the respective 
characteristic equations 



do 



P"(cos 9) 



= 



(3) 

(4) 



p;(cos e ) = o 

for a specified horn half angle 0„. 

These eigenvalues can be computed by a variety of numerical meth- 
ods. For example, one can represent the associated Legendre function 
in terms of the hypergeometric function' 1 as 



P:(cos d) = Cs\n m/2 d F (l + m + u, m - v, m + 1;- 9 °° s d ) 



(5) 
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where C is a constant. The v- or /i-zeros for a given 8 can then be 
found by a variety of root finding techniques. It is, however, worth 
noting that a first-order approximation may be determined from the 
formula 

P"(cos d) W J m (Vv(y + 1) 6) + 0(0 2 ), (6) 

where J m denotes a Bessel function. Since the roots of Bessel functions 
are well tabulated, 7 (6) can be conveniently used to provide either 
an estimate of the eigenvalue or a starting value for an iterative algo- 
rithm. The phase errors associated with this approximation will be 
discussed in a later section. 

To indicate the behavior of the zeros and to provide helpful in- 
formation for design of dual mode conical horns, a partial list of v 
and /x values, computed using (5), has been prepared and is given in 
the Appendix. 

III. BEHAVIOR OF THE RADIAL FUNCTIONS. PRECISE CALCULATION OF PHASE 
PROGRESSION. 

Having obtained the appropriate eigenvalues as defined by (3) and 
(4), we may then proceed to the more interesting calculation of the 
radial dependence. In principle, the phase shift between the two 
spherical surfaces r = n and r = r 2 is given by 

8a TE = arg [#%«] " arg [H^krJ] (7) 

-^•g{^[VxHl\\(x)]Y^ (8) 

for the TE and TM modes, respectively, where arg ( ) denotes the 
phase angle associated with a complex number. 

The difficulty which arises when one attempts to utilize these ex- 
pressions is essentially one of computation, due to the particular 
regime of order and argument frequently encountered in analysis of 
conical horn waveguides. We are particularly concerned here with the 
so-called transition region where the argument and order of the Hankel 
functions are large and comparable. For example, when O = 5°, v = 
20.6155 and /* = 43.4109 (see Appendix), which means that we must 
allow for a range of arguments increasing from these values. 

Asymptotic formulas for Bessel functions have, of course, been 
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studied in great detail. Asymptotic forms in the transition region kr 
J> v have been given, for example, by Watson and Langer. 8 Although 
of mathematical interest, these typically give only the limiting be- 
havior of the function, with the remainder specified within some order. 
However, engineering design generally requires more precise results, 
which can be obtained only with the aid of complete asymptotic ex- 
pansions. 

Expansions particularly appropriate for our problem have been 
given by Olvcr. 4 These formulas represent the Bessel functions as 
asymptotic series in terms of reciprocal powers of the order, and are 
valid asymptotic expansions for all values of argument. Although 
their derivation is very complicated, and will not be discussed here, 
we shall state the general form of the result and indicate several sim- 
plifications which are valid for most problems involving phase pro- 
gression in conical horns. 

The complete asymptotic expansion for the Hankel function of the 
first kind, following Olver's notation, is given by 



Hl"M 



4r y /AjfeJfl - iBijvh) 
k » 2 - + " ^ h~^ r y (9) 

In this expression, Ai, Ai', Bi, and Bi' represent Airy functions and 
their derivatives, 9 and £ is a constant related to x by the formula 

f = -!l(.f 2 - l) 5 - | sec" 1 *}*. (10) 

The coefficients A„{£) and /?„(£) are determined through an auxi- 
liary sequence {U„(t)} defined by the recursion formula 

U0(t) = \ , , W 

UJ® = r(1 2 n U' n -M + | jf (1 - 5t 2 )U n . 1 (t) dt, 

where the prime denotes differentation. The A n and B n can then be 
found using the relations 

Mt) = t ^iS^ (12) 

m-0 S 

s-(f) = -Z a :£ 2 :r; /2 +1 (t) , (is) 

m-0 i 
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whore r - (l-.r 2 )-* and {a,,,}, {b m } arc given by 

a = b a = 1 (14) 

, . m (2m + l)(2m + 3) ■ ■ • (6m - 1) ( . 

a m (m > 0) = , ,., nm U°J 

mV ' m! (144) 

6 m (wi > 0) = -7 — --7 a m . (16) 

6m — 1 

Since for the range x > 1, which is of interest here, £ is negative and 
t is imaginary, it is also necessary to define the proper branches, which 
are as follows: 

r = i(x 2 - 1) H (17) 



f* FF 



(18) 



(19) 



f 3m/2 " (-f) 3 "" 2 

Using a table of Airy Integrals, one can proceed to evaluate the 
required Hankel functions to whatever accuracy is needed. As an in- 
dication of the number of terms required in a typical calculation, it 
has been observed by J. A. Cochran and C. M. Nagel* that for v ^ 10, 
four decimal place accuracy can be obtained simply by using terms 
including B and A } . The coefficients required for most horn calcula- 
tions are thus given by 

A = 1 (20) 

. 81t 2 + 462rj + 385t? _ 7(3r t + 5r 3 ) 455 . . 

Al = ~ 1152 1152f* "•" 4608tf ^ ' 

_ 3ti + 5r, _ 5 / 22 ) 

Bo ~ 24f5 ml ' 

where Tl = (x 2 - 1) _i and & = -£. 

IV. APPLICATIONS 

In this section we shall discuss several applications of the preceding 
results. After presenting examples of phase progression for different 
modes we consider the phase errors introduced by approximating the 
eigenvalues /x and v. Finally, we examine what might be called the 

* Private communication. 
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"quasi-cylindrical approximation," in which phase progression is 
calculated by considering the horn to be a cylindrical waveguide with 
slowly varying cross section. 

4.1 Phase progression of the TEf t and TMf x modes. Effect of errors 
in v and n. 

A qualitative understanding of the phase progression properties of 
conical waveguide modes can be achieved by regarding a horn as a 
cylindrical waveguide with gradually increasing cross section. Although 
such a model has limitations, which will be discussed later in Section 
4.2, it correctly predicts the fact that the phase progression rates 
for both classes of modes begin at relatively low values and increase 
monotonically toward that corresponding to the far field of a spherical 
wave in free space. This limiting behavior is reached when the conditions 
kr » v or kr » /z are satisfied, corresponding to conical TE or conical 
TM modes, respectively. 

Fig. 2 shows, for example, a direct computation of the phase shift as 
a function of kr 2 — krj for the conical TE U mode for half angles O = 
3°, 10°. The value kr t is in each case taken to be that corresponding 
to the cut-off cross section of a cylindrical waveguide, i.e., 



kri = 



(23) 



sin O ' 
where £ ^ 1.84118. Fig. 3 shows analogous results for the TM< mode, 
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Fig. 2 — Phase shift for the TEf, mode relative to the cross section at which cut 
off would occur for a cylindrical waveguide. f TE » 1.841 IS. 
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Fig. 3 — Phase shift for the TM< mode relative to the cross section at which cut 
off would occur for a cylindrical waveguide. f TM ~ 3.83171. 

the principal difference being in the more gradual increase in phase 

near kr x . 

In view of the difficulty in computing the v- and /t-zeroes of the 
Legendre functions, as required by (3) and (4), it is of practical in- 
terest to determine how an error in the eigenvalue will effect the cal- 
culation of phase shift. If, for example, we denote by A TE the error in 
-phase shift due to a small error 8v in the eigenvalue, then from (7) we 



Table I — Maximum Phase Error Due to 0.10 Percent 
Miscalculation of v or /x 

k(r 2 - r,) = 100 
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ATE = (,./2) Sv 


ATE actual 


3° 

5° 

10° 


0.035 
0.021 
0.010 


3.15° 
1.89° 
0.90° 


2.04° 
1.26° 
0.61° 



TMii mode 
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«M 


atm =, (,/2) a„ 


ATM actual 


3° 

5° 

10° 


0.073 
0.043 
0.021 


6.87° 
3.87° 
1.89° 


4.68° 
3.14° 
1.16° 
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have that 

A TE = jl arg [Hl\\(kr 2 )] - £ arg [ffJV^i)]} p ' 5 "> ( 24 > 

where v denotes the correct eigenvalue. Although the calculation 
required by (24) is, in general, very difficult, it is relatively simple to 
obtain an upper bound to A TE . First, it can be shown (see Ref. 7, p. 
368), that at the cutoff radius r, , the argument kr x is approximately 
equal to the eigenvalue v (or fi) and furthermore that the partial 
derivative at that value tends toward zero as v or /x, becomes very 
large. It follows that an upper bound on the error A TE can be obtained 
by neglecting the second term on the right side of (24) and letting 
kr 2 — > oo. In this way we find, from the asymptotic behavior of the 
Hankel functions (see Ref. 6, p. 85) , that 

max | A TE | = I 8v | 5 (25) 

and, in a similar way for A™, 

max I A™ I - I fe 1 1. (26) 

In Table I we present a comparison between the actual computed 
error in differential phase shift, for an assumed relative 0.1 percent 
error in the eigenvalue, and the upper bound as determined by (25) 
and (26). The results indicate that the predicted estimates are quite 
reasonable. Note that the larger phase errors for smaller angles are due 
simply to the fact that the eigenvalue and hence the absolute error 
is greater. 

The principal purpose of the previous exercise was to determine 
what error might be expected from using (6), which expresses the 
Legendre function in terms of a Bessel function. Results for the TM U 
mode show that (6) is sufficiently accurate in predicting the /^-zeroes 
that the maximum differential phase error for horns up to 30° in half 
angle should be less than 1°. As might be expected, however the same 
approximation applied to the v-zeroes of the derivative of the Legendre 
function is not as accurate. Nevertheless, as shown in Fig. 4, for a horn 
with half angle equal to 30° the maximum phase error is approxi- 
mately 6°, which would ordinarily be acceptable. 

4.2 Evaluation of the Quasi-Cylindrical Approximation 

The difficulty of making precise calculations of phase progression in 
conical horns has led to the use of an approximate formula derived 
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Fig. 4 — Maximum error in phase shift using Bessel approximation to the Legendre 
function. 

by assuming that the horn behaves as a cylindrical waveguide with 
gradually increasing cross sectional radius. The phase shift is then 
determined simply by integrating the local waveguide propagation 
constant, with the result that 



tf 



- J — P V(fr 

tan O Lr, 



;?• sin d ) — r — f cos" 



kr sin 



;]■ 



(26) 



where f is the characteristic value for the particular mode (e.g., for 
the TE,° mode r ~ 1-841 and for the TM° mode f ^ 3.832). This 
formula is, in fact, asymptotic to the true phase shift in certain limits. 
For example, let (kr sin O ) be fixed and let O ~> 0. Then it can be shown 
that (26) becomes essentially equivalent to the simple Debye approxi- 
mation. (See Ref. 7, p. 366.) However, this formula is known to be 
invalid in the range where order and argument are comparable. Never- 
theless, it is useful to investigate the properties of (26) in order to 
determine what errors accompany its use. 

Fig. 5 shows the resulting error in differential phase shift when the 
quasi-cylindrical approximation is applied to the TE< mode, with 
kr x = f/sin O (corresponding to the cut-off diameter) and O = 5°. 
The error is seen to grow very rapidly at first, showing that the quasi- 
uniform approximation predicts too slow an increase in propagation 
constant with increasing cross section. Eventually, the error curve 
approaches a linear variation. This asymptote can actually be predicted 
fairly well by using the large argument behavior of the Hankel function, 
combined with the fact that for large v (see Ref. 7, p. 36S), 

argff^W^ -t/S. (27) 
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By letting kr sin $ — » °o in (26) and using (6), which relates £ to v, 
we find that in the limit of large kr 



8a TE - 5/3 TK 
which, for small 6„ reduces to 



12 fTK U, tan 6 J 2 



E+[l -fi^Jfc. 



tan O 



X TE sflTE _, T rTK^T . , 0<"> , 



(28) 



(29) 



This result, shown in Fig. o as the dotted line, is seen to predict very 
accurately the asymptotic behavior of the error. 

Fig. 6 shows the corresponding error for the TM,* mode. In con- 
tradistinction to the previous example, the quasi-cylindrical approxi- 
mation at first predicts too high a phase progression rate, but eventually 
also conforms to the linear error predicted by the last term of (28). 
The formula analogous to (29) is given by 



Stm* 1 



8a - 5/3 = -— -— d Q + — kr. 



(30) 



This result, shown in Fig. 6 as the dotted line, is also seen to cor- 
rectly describe the asymptote. 

A salient feature of these results is that the linear portion of these 
curves is quite independent of the type and order of the mode being 
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Fig. 5 — Error in phase shift due to quasi-cylindrical approximation for the 
TE< mode. f TB « 1.84118. 
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Fig. 6 — Error in phase shift due to quasi-cylindrical approximation for TMf, 
mode. Ttm ~ 3.83171. 



considered. This immediately implies that the error in differential 
phase shift between modes resulting from the quasi-cylindrical ap- 
proximation is always bounded at a value easily predicted by (29) 
and (30). This latter result is considered to be one of the more sig- 
nificant conclusions of this study. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have considered the phase progression properties of 
conical waveguide modes. The principal difficulties have been in com- 
puting Bessel functions over their so-called transition region. It is 
suggested that, in view of the typically large orders involved, the 
asymptotic expansions due to Olver are the most applicable. An exam- 
ination of the quasi-cylindrical approximation has shown that this 
latter formula, although not necessarily accurate for evaluating phase 
progression of a particular mode, can be used to determine differential 
phase shift between modes with an error which is bounded over the 
conical region. 
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APPENDIX 



Roots of PI (cos d) = and d/dd P\ (cos 6) = for a Given Value of 

In order to assist in design of dual-mode conical horns, we have 
prepared an accurate table of the roots of the associated Legendre 
Function and its derivative for a given value of the horn half angle 0. 
These values of v and n ma}' then be used to evaluate the vector wave 
functions which characterize the propagation of the TEJ\ and TMf x 
modes in the horn section. 

The computational program consisted of using the hypergeometric 
series representation for the Legendre function, and determining the 
zeroes by a standard root finding method. The program was termi- 
nated when the value of the function was less than 10~ n in amplitude. 



Table II — Roots of P* (cos d) = 



e 


^ 


9 


n 


3.00° 


72.6819 


17.00° 


12.4239 


3.50 


62.3379 


17.50 


12.0552 


4.00 


54.3874 


18.00 


1 1 . 7070 


4.50 


4S.2893 


18.50 


11.3777 


5.00 


43.4109 


19.00 


11.0657 


-).-)() 


39.4196 


19.50 


10.7697 


6.00 


36.0935 


20.00 


10.4885 


6.50 


33.2792 


20.50 


10.2211 


7.00 


30.8669 


21.00 


9.9664 


7.50 


28.7764 


21.50 


9.7235 


N.00 


26.9471 


22.00 


9.4918 


8.50 


25.3332 


22.50 


9.2703 


9.00 


23.8985 


23.00 


9.0585 


9.50 


22.6110 


23.50 


8.8557 


10.00 


21.4597 


24.00 


8.6613 


10.50 


20.4146 


24.50 


8.4749 


1 1 . 00 


19.4645 


25.00 


8.2960 


1 1 . 50 


18.5970 


25.50 


8.1241 


12.00 


17.8019 


26.00 


7.9589 


12.50 


17.0704 


26.50 


7.7998 


13.00 


16.3952 


27.00 


7 . 6467 


13.50 


15.7700 


27.50 


7.4992 


14.00 


15.1894 


28.00 


7.3570 


14.50 


14.6490 


28.50 


7.2197 


15.00 


14.1446 


29.00 


7.0871 


15.50 


13.6728 


29.50 


6.9591 


16.00 


13.2304 


30.00 


6.8354 


16.50 


12.8149 
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Table III- 


—Roots of 


d/dd P? 


(COS 6) = 


e 
3.00 


34.6743 


e 
17.00 


5.7637 


3.50 


29.6526 


17.50 


5.5881 


4.00 


25.8867 


18.00 


5.4224 


4.50 


22.9581 


18.50 


5.2657 


5.00 


20.6155 


19.00 


5.1174 


5.50 


18.6992 


19.50 


4.9768 


6.00 


17.1026 


20.00 


4.8432 


6.50 


15.7518 


20.50 


4.7163 


7.00 


14.5943 


21.00 


4.5955 


7.50 


13.5913 


21.50 


4.4804 


8.00 


12.7139 


22.00 


4.3706 


8.50 


11.9400 


22.50 


4.2658 


9.00 


11.2522 


23.00 


4.1656 


9.50 


10.6370 


23.50 


4.0697 


10.00 


10.0835 


24.00 


3.9779 


10.50 


9.5828 


24.50 


3.8900 


1 1 . 00 


9.1279 


25.00 


3.8056 


11.50 


8.7126 


25.50 


3.7246 


12.00 


8.3321 


26.00 


3.6467 


12.50 


7.9822 


26.50 


3.5719 


13.00 


7.6593 


27.00 


3.4999 


13.50 


7.3605 


27.50 


3.4306 


14.00 


7.0831 


28.00 


3.3638 


14.50 


6.8250 


28.50 


3.2995 


15.00 


6.5842 


29.00 


3.2374 


15.50 


6.3591 


29.50 


3.1775 


16.00 


6.1481 


30.00 


3.1196 


16.50 


5.9500 
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